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tron  Act  of  19^9 . 
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ABSTRACT 


The  loss  for  a  partially  magnetized  material  can  be  characterized 
by  n"  of  the  completely  demagnetized  state  if  u^/to  —  0.  7.  Results  are  pre¬ 
sented  for  a  number  of  yttrium  iron  garnets  with  aluminum  and  gadolinium 
substitutions  and  magnesium  manganese  spinels  with  aluminum  substitutions. 
|i"  has  been  measured  at  frequencies  of  5.  5,  7.3,  9.0,  11.4,  13.8,  16.9GHz. 
The  results  can  be  approximated  by  a  relation  jr1  -  A(u^/u)  .  The  YA1G 
materials  are  low  in  loss,  MgMnAlF  slightly  higher,  and  the  YGDA11G  con¬ 
siderably  higher  in  loss.  Two  theoretical  modes  have  been  developed  which 
are  an  improvement  over  Rado's  theory  and  give  qualitative  descriptions  of 
the  partially  magnetized  permeability  tensor. 

High-power  parallel- pump  measurements  have  been  made  on  the 

same  materials  at  frequencies  of  5.  5  and  9.  2  GHz  and  at  fields  of  0,  4ttM  /3, 

s 

and  H  .  .  These  materials  have  the  same  ranking  with  respect  to  threshold 
min  & 

as  they  have  for  low- power  loss. 
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SECTION  I 


INTRODUCTION 


Modern  radar  and  communication  systems  are  putting  requirements 
on  microwave  ferrite  devices  which  are  a  test  of  the  ingenuity  of  the  device 
designer  and  the  materials  physicist.  Phased- array  systems  require  tem¬ 
perature  stability,  cheapness,  reproducibility,  and  light  weight.  Low  inser¬ 
tion  loss  and  high  peak- power  performance  are  highly  desirable.  The 
driving  circuits  must  be  simple,  dependable,  and  inexpensive.  There  is  also 
a  trend  to  miniaturization  and  integration.  In  any  given  application,  it  is  im¬ 
possible  to  maximize  every  single  one  of  these  performance  characteristics 
simultaneously  and  therefore  some  sort  of  a  trade-off  procedure  is  required. 

In  order  to  be  able  to  make  a  performance  trade-off,  it  is  necessary 
to  be  able  to  predict  the  effect  of  design  procedures  upon  performance  char¬ 
acteristics.  The  goals  of  the  present  program  are  to  study  how  the  perme¬ 
ability  of  partially  magnetized  materials  depends  on  material  parameters 
and  to  catalogue  the  loss  properties  of  a  number  of  the  important  phase 
shifting  materials.  It  should  then  be  possible  to  use  computer-aided  design 
procedures  in  the  trade-off  procedure. 

The  results  of  the  last  half  year  are  reported  in  two  sections.  The 
first  section  deals  with  th^  low  power  loss  parameter  p"  of  the  demagnetized 
state  as  measured  on  a  number  of  materials  across  a  range  of  frequencies. 
Previous  work  has  determined  the  dependence  of  real  parts  of  the  perme¬ 
ability  tensor  upon  material  parameters.  This  permits  phase-shifiting  calcu¬ 
lations.  The  present  work  permits  the  calculation  of  insertion  loss.  The  sec¬ 
ond  section  deals  with  the  high-power  performance  of  these  same  materials. 


SECTION  II 


PERMEABILITY  TENSOR  OF  PARTIALLY  MAGNETIZED 
MAGNETIC  MATERIALS 


One  of  the  purposes  of  this  contract  has  been  to  describe  the  Joss 
part  of  permeability  tensor  for  a  wide  variety  of  materials  in  partially  mag¬ 
netized  states.  Earlier  measurements  of  p11  and  k"  about  the  hysteresis 
loop*  have  shown  that  when  <  0.  7,  «"  «  p".  This  is  confirmed  theo¬ 

retically  by  calculations  of  Schlomann  described  in  an  appendix  to  this  report. 
In  addition,  it  was  found  that  p"  is  relatively  insensitive  to  position  on  the 
hysteresis  loop  in  this  region.  Only  a  shallow  minimum  at  4rrM  =  0  is  ob¬ 
served.  Measurements  of  p"  have  shown*  that  it  has  a  maximum  at  4ttM  -  0 

z 

where  it  equals  p".  It  is  thus  reasonable  to  characterize  the  partially  mag¬ 
netized  states  by  a  single  loss  constant  (when  u^/u  <  0.7),  namely  p:l  at 
4itM  =  0  which  is  nearly  an  upper  limit  on  all  of  the  loss  terms. 


Measurements  of  p'1  in  the  demagnetized  state  have  been  made  on 
seven  garnets  with  various  amounts  of  aluminum  and  gadolinium  substitution 
and  on  five  magnesium  ferrites  with  various  amounts  of  aluminum  substitu¬ 
tion.  The  samples  had  magnetizations  ranging  from  725  to  2475  gauss. 

Most  of  the  samples  were  measured  at  six  frequencies  from  5.  5  to  16.7  GHz. 
This  covered  u^/u  ratios  from  0.  11  to  0.  89.  The  measurements  were  made 
on  thick  slabs  in  a  rectangular  cavity  as  described  in  Ref.  1.  Thick  slabs 

rather  than  the  more  conventional  rod-  or  sphere-shaped  samples  were  used 

-4 

to  improve  sensitivity  in  the  region  p' '  <  10  .  The  results  are  displayed 

in  Fig.  1  on  a  log-log  plot  of  p"  vs  u^/u.  The  experimental  accuracy  is 
believed  to  be  comparable  to  the  scatter  from  the  interpolated  curves, 
namely  ±  (10  percent  +  0.  00001).  The  data  follows  a  simple  power  law 
relationship: 


p"  =  a  (!M)N 

rdemag  u 


(1) 


for  all  u^/u  <  O'  7.  At  greater  ratios  of  co^/ to,  p"  increases  more  rapidly, 
becoming  extremely  large  as  u^/to  approaches  1.  0.  The  measured  values 
of  A  and  N  are  listed  in  Table  I.  It  should  be  noted  that  increases  in  gado¬ 
linium  content  in  YIG  and  increases  in  aluminum  content  in  MgMn  ferrite 
increase  A  and  decrease  N,  whereas  an  increase  in  aluminum  content  in  YIG 
has  little  effect  on  loss  other  than  through  the  change  in  magnetization. 

For  a  saturated  sample  at  zero  field 


-  0 


vAHeff  UM 

- 2  Z 

u 


(2) 


where  AH^  =  2/yT  is  an  effective  linewidth  proportional  to  the  intrinsic 

relaxation  rate  of  the  material,  not  the  ferromagnetic  resonance  linewidth. 

A  subthreshold  parallel  pump  loss  at  zero  field  has  been  derived  by  Joseph 
2 

and  Schldmann.  Their  results  show  that  |jY  ^  q  depends  on  tempera¬ 

ture,  and  the  exchange  constant,  but  not  on  AHg^.  In  the  demagnetized 
state,  1/3  of  the  domains  are  being  parallel  pumped  and  2/3  of  the  domains 
are  being  perpendicularly  pumped.  We  may  thus  assume  that  p"  in  the  de¬ 
magnetized  state  with  u„/u  well  below  1. 0  is  given  by: 

p"  \  p"  n+U"  „  n  •  (3) 

rdemag  3  H=0  3  z,H~0 

The  term  l/3  p"  n  is  plotted  in  Fig.  1  for  4ttM  =  800  and  1  7  50  gauss 
z,  H-U  _q  s  2 

using  the  exchange  constant  of  YIG  (D  =  5  •  10  Oe-cm  ). 

From  our  results  it  appears  that  for  the  very  low-loss  materials 
G-  113,  G-810,  andG1210,  i.  e.  ,  those  having  the  lowesi  value  of  A,  a  sig¬ 
nificant  part  of  HL^ernag  arises  from  the  subthreshold  parallel  pump  loss.  For 

the  higher  loss  materials  the  p-,1  arises  almost  entirely  from  the  usual 

6  rdemag  J 

perpendicular  pump  loss. 
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TABLE  I 


LOSS  PARAMETERS  OF  GARNETS  AND  MAGNESIUM  FERRITES 


Material 

Composition 

4irM 

S 

(gauss) 

A 

N 

AH  ,,  at 
eff 

uM/w=1-° 

(Oe) 

G- 113 

YIG 

17  50 

-4 

7.  5  •  10 

2.  5 

3.0 

G-  1210 

YIG  +  A1 

1200 

10.  0 

2.4 

3.0 

G-  1200 

YIG  +  A1  +  Gd 

1200 

34.  0 

1.9 

12.  0 

G-  1001 

YIG  +  Gd 

1200 

50.  0 

1.  8 

18.  0 

G-810 

YIG  +  A1 

800 

8.  0 

2.  8 

1.5 

G-800 

YIG  +  A1  +  Gd 

800 

43.  0 

1.9 

15.  0 

G-  1 005 

YIG  +  Gd 

725 

95.  0 

1.3 

21.0 

MMF-10 

MgMn  ferrite 

2475 

20.  0 

2.  0 

15.  0 

1-390 

MgMn  ferrite 
with  second  phase 

2200 

16.  0 

2.2 

11.0 

1-105 

MgMn  ferrite  +  A1 

1750 

18.  0 

2.  0 

9.0 

1-109 

MgMn  ferrite  +  A1 

1300 

28.  0 

1.  9 

11.0 

1-414 

MgMn  ferrite  +  A1 

750 

45.  0 

1.7 

10.  0 

Combining  Eqs.  (1-3)  we  may  set 


2  vAHeff  WM  ,  1 

7  - 2 - r  +  7  ^z,  H=0 

u 


to 


A(- 


_M 

to 


N 

) 


or 


AH  ,,  =  4irM 
eff  s 


(— ) 
to 


!  N 

>  3A  (— ) 

\  to 


'  h 


z,  H=0 


(4) 


(5) 


When  the  parallel  pump  term  is  small,  (N-2)  gives  the  frequency  dependence 

of  AH  The  values  of  AH  ..  at  to.  ./to  -•  1.0  are  easily  obtained  from  Eq.  (5) 
eff  eff  M'  J  n 

and  are  tabulated  in  Table  I. 


Measurements  of  p'  on  the  same  materials  were  made  at  each  fre¬ 
quency  at  the  same  time  as  the  p"  measurements.  The  measured  values  of 
p1  are  plotted  vs  tO-^/to  in  Fig.  2,  and  it  can  be  seen  that  the  points  for  all 
the  materials  fall  close  to  a  single  curve.  Also  included  in  this  figure  are 
some  of  the  measurements  of  p'  on  demagnetized  magnesium  ferrite  at  a 
single  frequency,  5.  5  GHz,  reported  in  Kef.  1.  In  these  measurements, 
temperature  (and  hence  ioM)  was  varied. 


A  theory  of  the  partially  magnetized  state  proposed  by  Sandy  is 

3 

appended  to  this  report.  It  extends  Rado's  earlier  theory  of  the  partially 
magnetized  state  to  include  correlation  between  internal  demagnetizing 
fields  and  domain  orientations.  This  results  in  substantially  improved 
relations  for  p'  and  p'  vs  -y4ttM/w  and  y4ttM  /u.  In  particular  for  the 
demagnetized  state  it  yields 


=  1- 


\  N  (ujm  +  YHa/N)2/w 
i  -  An2  (0JM+yHa/N)2/uj2 


(6) 


where  N  is  an  average  demagnetizing  factor  and  H  is  an  average  aniso- 

cl 

tropy  field.  The  best  fit  to  the  data  in  Fig.  2  for  Ujyj/w  <  0,  9  is  obtained 


vs  at  4ttM  =  0  for  Various  Garnets  and  Mg  Ferrites 


SECTION  III 


THE  HIGH  POWER  PROPERTIES  OF  PARTIALLY 
MAGNETIZED  MATERIALS 

The  same  materials  whose  properties  have  been  reported  in  the 

previous  permeability  section  are  also  being  measured  at  high  power  levels. 

The  program  calls  for  parallel  pump  measurements  at  frequencies  of  3.0, 

5.  5,  9.2,  and  16.  5  GHz  at  dc  fields  of  0,  4ttM  /3,  and  H  .  ,  where  H  . 

s  min  min  4 

is  the  dc  field  for  the  minimum  parallel  pump  threshold  (butterfly  minimum). 
Spherical  samples  and  cavity  perturbation  techniques  are  used.  For  spheri¬ 
cal  samples  H^c  =  0  corresponds  to  complete  demagnetization  and  H^c  = 

4ttM  /3  corresponds  to  complete  magnetization  with  an  internal  dc  field  ap- 
s 

proximately  equal  to  zero.  Measurements  have  been  completed  at  5.  5  and 
9.  2  GHz.  The  results  are  summarized  in  Tables  II  and  III.  The  tables  in¬ 
clude  the  spinwave  linewidth  AH^  computed  for  H^  =  4itM/3  and  using 

the  parallel  pump  threshold  equation4 


h  .  =  — 
crit 


where  u  is  the  signal  frequency  and  =  y4TTMg  with  y  =  2tt  times  the  gyro- 
magnetic  ratio.  In  order  to  compare  the  strength  of  the  relaxation  mech¬ 
anism  responsible  for  low  power  insertion  loss  with  that  responsible  for  the 
high  power  threshold  the  effective  linewidth,  AH^  described  in  the  pre¬ 
vious  section,  is  included  in  the  table.  AH  ^  was  computed  for  the  same 
signal  frequency  as  the  parallel  pump  threshold  experiment  was  performed. 

The  following  observations  can  be  made  on  the  data  presented  in 
Tables  II  and  III.  The  YIG  (G1 13)  and  the  YA1IG(G810,  G1210)  have  the  low¬ 
est  AH^  and  AH^  among  the  garnets.  The  addition  of  Gd,  YGDA1IG  (G800, 
G1200),  YGdIG  (G1005,  G1 001)  causes  a  significant  increase  of  AH^  and 
AH  .  The  MgMnFe  (TT1-390)  and  the  MgMnAlFe  (TT1-414,  TT1-109,  and 
TT1-I05)  are  all  quite  similar  in  their  AHj  and  AHg^.  For  those  materials 
which  contain  Gd,  AH  „  is  as  much  as  factor  of  two  larger  than  AH,  .  This 
is  due  to  the  frequency  dependence  of  the  rare-earth  relaxation  mechanism 


9 


(i.e.  ,  the  spinwave  is  excited  at  w/2,  and  the  uniform  precession  is  excited 
at  u).  The  YIG  and  YA1IG  are  the  lowest  loss  materials  available.  The  fact 
that  this  particular  specimen  of  G810  had  the  lowest  threshold,  lower  than 
that  of  the  G 1 1 3 ,  is  an  example  of  microstructural  effects  (e.  g.  ,  values  of 
AH^.  as  low  as  0.7  Oe  have  been  seen  on  some  Raytheon  YIG).  Therefore  it 
must  be  emphasized  that  the  values  of  AH^  and  AH  ^  appearing  in  these 
tables  are  neither  upper  or  lower  limits,  but  rather  typical  values  for  ma¬ 
terials  obtained  from  a  commercial  supplier.  A  change  of  microstructure 
by  an  alteration  of  the  material  processing  procedure  will  give  different 

values  of  AH.  and  AH  The  difference  in  AH,  for  H  ,  =  4ttM/3  and  H 

k  eff  k  dc  min 

is  due  to  a  wavelength  dependence  in  AH^* 

A  desirable  feature  is  low  insertion  loss  and  high  peak  power  thresh¬ 
old  or,  in  terms  of  the  relaxation  mechanisms,  AH,  »  AH  ...  In  the  case 
of  the  YA1IG  at  least,  AH^  >  AH^.  For  the  gadolinium  materials  and  the 
spinels  AHeff  >  AH^.  From  rare-earth  relaxation  theory  one  expects 
AHeff(u)  =  2 AH^(u/2)  because  of  the  frequency  dependence.  Since  no  theory 
exists  which  accounts  for  the  microstructural  effects  the  only  comments 
which  can  be  made  are  that  the  microstructure  these  particular  YA1IG 
samples  has  a  stronger  effect  on  AH^  than  on  AH^,  while  on  these  particu¬ 
lar  MgMnAlFe  samples  the  opposite  is  true.  The  relation  between  AH^  and 
AHgff  f°r  those  samples  containing  gadolinium  is  what  would  be  expected 
from  rare-earth  relaxation  theory. 
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APPENDIX  1 


PERMEABILITY  OF  PARTIALLY  MAGNETIZED  FERRITES 


Frank  Sandy 
Raytheon  Company 
Research  Division 
Waltham,  Massachusetts 


ABSTRACT 


Existing  theories  of  the  partially  magnetized  state  incorrectly 

predict  that  pL*  =  (jl1  =1.0.  This  arises  from  their  neglect  of  the  spatially 
z 

fluctuating  internal  rf  magnetic  field  produced  by  rf  magnetic  poles  on 
domain  and  grain  boundaries.  This  report  attempts  to  take  these  fields 
into  account  in  a  statistical  manner.  Although  the  model  used  is  quite 
crude  the  theory  gives  good  agreement  for  [P  and  vs  um/w  in  the  de¬ 
magnetized  state,  qualitatively  predicts  their  dependence  on  4ttM,  and 
correctly  predicts  hysteresis  in  permeability  vs  4ttM.  It  fails  however 
to  yield  a  significantly  better  approximation  to  k'  vs  4ttM  than  the  earlier 
theories. 


i 


I.  INTRODUCTION 


Existing  theories  1  of  the  real  part  of  the  permeability  of  par¬ 
tially  magnetized  ferrites  give 

H'  =  1 

H-1  =  1 
z 

k'  =  y4TrM/id  .  (1) 

The  equation  of  k'  is  qualitatively  correct  but  the  equations  for  h1'  and  jjl' 

z 

are  not,  because  the  derivations  do  not  properly  take  into  account  the  inter¬ 
nal  rf  magnetic  fields  produced  by  rf  poles  on  domain  boundaries.  These 
fields  are  spatially  fluctuating  and  average  to  zero  over  the  whole  sample 
but  still  contribute  to  the  permeability  of  the  material.  When  >  1. 

Polder  and  Smit3  showed  that  the  effect  of  these  fields  was  far  from  negli¬ 
gible  and  was  in  fact  the  source  of  what  is  known  as  low  field  loss.  In  the 
region  ^  1.  the  effect  of  these  fields  is  not  as  spectacular  but  is  still 

responsible  for  the  errors  in  Eq.  (1). 

A  proper  statistical  analysis  of  these  spatially  fluctuating  internal 
fields  is  not  available  at  this  time.  It  is  the  purpose  of  this  article  to  show 
by  means  of  an  oversimplified  model  how  these  fields  affect  the  permeability. 
The  numeric  results  however  do  not  fit  the  data  perfectly  indicating  the 
necessity  to  develop  a  more  sophisticated  model  of  a  partially  magnetized 
ferrite.  However,  the  do  show  most  of  the  qualitative  features  of  the  de¬ 
pendence  of  p.1 ,  |jl'  ,  and  k'  on  y4ttM/cj  and  y4irM  /u. 

z  s 

H.  NOTATION 

Before  describing  the  physical  model,  it  is  necessary  to  discuss 
the  notation  that  will  be  used  to  describe  the  various  components  of  mag¬ 
netization  and  magnetic  field.  The  notation  and  approach  are  similar  to  that 
of  Rado.  We  let  H,p  be  the  total  internal  magnetic  field  at  any  point  in  the 
sample. 
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where  H  is  the  dc  component  of  Hrp  and  E  is  its  rf  component.  H  consists 
of  an  applied  magnetic  field  less  dc  demagnetizing  fields  from  the  surface, 
anisotropy  fields,  and  demagnetizing  fields  from  dc  poles  on  domain  bound¬ 
aries.  It  is  often  assumed  that 

H  «  4ttM  (3) 

s 

where  4ttMs  is  the  saturation  induction  of  the  material.  The  total  mag¬ 
netization  M„  of  the  sample  is  composed  of  a  dc  part  M  and  a  small  rf 
X  s 

component  m.  We  can  write 

Mrp  =  M  +  m  (4) 

T  s 

where  the  magnitudes  of  M  and  M  are  equal  since  we  are  dealing  only 

S  1 

with  the  low  power  limit.  Each  of  the  quantities  H,  M  ,  and  E  can  then 
be  divided  into  an  average  part  denoted  by  brackets  ^  )  and  a  spatially 
fluctuating  part  denoted  by  6. 

H  =  (h)  +  6H 
h  =  ^h^  +  <5h 
m  =  +  6m 

MB  =  <Ms>  +  6M  .  (5) 

We  shall  also  denote  the  average  or  measured  magnetization  by  its 

usual  symbol,  M.  The  average  internal  magnetic  field  ^h^>  is  just  the 
applied  rf  magnetic  field  less  the  rf  demagnetizing  field  from  the  sample 
boundary.  The  permeability  tensor  p 


which  we  are  trying  to  find,  is  defined  by 


4tt  (n\^  =  (|i  -  1 )  •  ^h/1 


III.  THEORY 


We  shall  obtain  [i  by  starting  from  the  torque  equation. 


where 


Mt  =  jum  =  -y(MT  X  Ht) 


=  g  le  1 


By  incorporating  Eqs.  (2),  (4),  and  (5)  into  (8)  making  use  of  the  approxi¬ 
mation  in  Eq.  (3)  and  neglecting  terms  quadric  in  the  rf  fields,  we  obtain 

(m)  +  6m=  -X  {(Ms)  X  (K)  +  (  Mg)  X  SK  +  5 M  x(E)  +  6 M  X  6E  }  .  (JO) 

Taking  the  spatial  average  of  Eq.  (10)  and  noting  that  by  definition 

(6m)  =  (6h)  =  (6m)  =  0  (11) 


we  get 


W  =  jjf  (  (ms>  X  (k)  +  (sm  X  al>  } 
=  ^  {  M  X  (h)  +  ^6M  X  6E  )} 


At  this  point  Rado  asserts  that  6M  and  6E  are  uncorrelated  so  that 


(6M  X  6h)  =  (6M>  X(6E)=  0 
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and  hence 


m  =  -  X  M  X  (E)  .  (14) 

If  we  let  M  be  in  the  z  direction  and  let  ^h^be  in  the  x,  z  plane  (where 

/*  s\ 

x  and  z  are  unit  vectors),  then 
m  =0 

m  <h\ 

y  J  w  \  x/ 

m  =  0  (15) 

z 

From  this,  one  obtains  Eqs.  (1).  Unfortunately  6M  and  6E  are  correlated 
and  Eq.  (13)  is  not  valid.  To  see  why  6M  and  6E  are  correlated,  we  must 
consider  a  model  that  gives  us  more  information  about  the  relationship  of 
6M  to  6E. 

IV.  MODEL 


We  shall  assume  that  each  domain  is  ellipsoidal  with  a  demag¬ 
netizing  tensor  N.  We  will  then  have  a  contribution  to  6h  from  the  rf  poles 
on  its  own  domain  boundary.  We  will  denote  this  by  6h^: 

6hj  =  -  N  •  (m  -  mj)  =  -  N  •  (m  -  ((m)  +  6m ^)) 

=  -  N  '  6m  +  N  -  6m  ^  (16) 

where  m^  is  the  rf  magnetization  of  the  material  surrounding  the  given 
domain.  We  shall  assume  that  6m^,  the  deviation  of  m^  from  ^  is 
uncorrelated  with  the  6M  of  the  domain  in  question.  We  will  also  assume 
that  the  contributions  to  6h  from  poles  on  all  other  domain  walls  other  than 
that  surrounding  the  domain  in  question  are  uncorrelated  with  6M.  Thus 
we  will  write 
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6K  =  -  N  •  <5m  +  terms  uncorrelated  with  6m.  (17) 

Since  we  will  be  concerned  only  with  (6M  X  6h^  ,  we  shall  neglect  all 
terms  uncorrelated  with  6M  at  this  time.  Thus  we  have 

fh  =  -  N  '  6m  .  (18) 

We  further  assume  without  any  justification  that  the  deviation  of  N  from 
its  average  is  uncorrelated  with  6M.  This  approximation  is  the  primary 
reason  for  what  disagreement  there  will  be  between  theory  and  experiment. 
However,  these  approximations  are  necessary  to  make  the  mathematics 
tractable.  We  thus  have  as  the  approximation  to  be  used  in  our  model 

6h  ~  -  4iTN6m  (19) 

where  N  is  an  average  scalar  value  of  the  demagnetizing  factor.  Although 
this  has  the  appearance  of  asserting  that  6h  is  related  to  6m  by  a  point 
relationship,  which  would  be  quite  dubious,  it  is  only  being  assumed  that 
that  part  of  5h  which  is  correlated  to  6M  is  related  to  6m  by  such  a  point 
relationship. 

Before  evaluating  ^6M  X  6h^  explicitly,  it  is  worthwhile  to  demon¬ 
strate  qualitatively  for  the  simple  case  where  4irM  =  0  just  why  (SM  X  6h) 
is  not  zero.  Rado's  argument  that  it  is  zero  states  that  for  any  domain  with 
a  given  6M  and  6h  there  is  another  domain  with  equal  probability  having  the 
same  6h  and  opposite  6M.  But  let  us  see  what  actually  happens  when  we 
reverse  6M  of  one  domain  keeping  everything  else  constant.  Since  4ttM  =  0, 
each  6M  has  the  magnitude  of  M  and  the  rf  magnetization  of  each  domain 
precesses  about  the  direction  of  6M  If  we  reverse  the  direction  of  6lVl, 
we  reverse  the  direction  of  precession  of  the  magnetization  and  also  the  di¬ 
rection  of  precession  of  the  6h  induced  by  the  magnetization.  The  components 
of  6h  parallel  to  the  drive  field  are  equal  for  these  two  oppositely  magnetized 
domains  and  the  components  of  6H  perpendicular  to  the  drive  field  and  to  6M 


will  be  opposite  in  sign.  Thus  the  components  of  <5M  X  6E  parallel  to  the 
applied  field  will  add  and  only  the  components  perpendicular  to  the  drive 
field  and  6M  will  cancel  each  other.  Thus  ^6M  X  <5lT)  has  a  component 
parallel  to  the  drive  field  yielding  a  contribution  to  p1  and  p1  not  predicted 
by  Rado. 

V.  EVALUATION  OF  (SM  X  6h) 

We  can  now  proceed  to  evaluate  ^6M  X  6h)  in  terms  of  our  model. 

We  first  rewrite  Eq.  (10)  including  the  terms  in  ^R\  and  <5H  which  were 

neglected  but  noting  that  M  X  H  =  0. 

s 

(m)  +dm  =  4{IX  ^ 'E )  +  M  X  6h  +  6M  X  (h)  +  6M  X  <5h 


+  (m)  X  +  (m)  X  6H  +  6m  X  (h)  +  6m  X  6H } 


Taking  the  spatial  average  of  both  side  as  before 


(20) 


(m)  =-X{Mx(H)-  (h)  X  (m)  +  ^6M  X  6E)  -  (6H  X  6m)  }  .  (21) 

J  ^  / 

Making  use  of  our  basic  assumption,  Eq.  (19),  we  can  rewrite  this  as 

(m)  =  -X  {M  X  (e)  -  (h)  X<m)  +  ((  6M  +  ^)  X  6h)  .  (22) 


Equation  (22)  has  two  additional  terms  that  were  absent  from  Eq.  (12).  The 
term  in  ^h)  X  ^m)  states  that  when  the  material  approaches  saturation 
and  the  spatially  fluctuating  terms  6M  and  6h  go  to  zero,  the  sample  be¬ 
haves  as  if  it  were  in  a  dc  field  of  ^h)  .  In  this  case,  however,  <^h) 
is  not  primarily  an  applied  dc  field  but  rather  an  average  of  anisotropy  and 


other  internal  fields.  For  |y  ^H)  |  =  «  w  this  term  contributes  an 

amount  -  to  p1  but  makes  no  contribution  to  k1  or  (J.V  The  next 


additional  term  6H/N  X  6h  has  an  effect  similar  to  6M  X  6h  only  smaller. 
Since  it  greatly  complicates  the  calculations  without  making  a  substantial 
contribution,  we  shall  calculate  (dM  X  6h^  assuming  H  =  0  and  replace  6M 
in  the  final  result  by  <5M  +  <5H/N.  With  H  =  0  subtracting  Eq.  (21)  from 
Eq.  (20)  yields 

6m  =  -  X  { M  X  6h  +  6M  X  (h)  +  6M  X  6h  -  (fiM  X  ah')  (23) 
Combining  the  first  and  third  terms  of  Eq.  (23)  and  makinguse  of  Eq.  (19)  we  get 


6h  =  {Mg  X  6Ti  +  6M  X  (h)  -  (SM  X  6H) } 


(24) 


We  must  now  solve  this  equation  for  (6M  X  6h^  in  terms  of 
4TT-yMs/u,  y  ^h^) ,  N,  and  the  statistics  of  6M.  Having  done  this, 

the  result  can  be  substituted  into  Eq.  (22)  and  the  components  of  (in')  found 
as  a  function  of  the  components  of  ^h^  .  The  details  of  this  algebra  and  the 
mathematical  approximations  necessary  are  described  in  the  appendix-  The 
results  of  this  calculations  are 


n’  - 


-  (ig')2  <61^))  _ 

{1  -  (ip.j  N2(M2  +  (5M2))}  w 


(25a) 


-  (i!Y>  N(  (6M2V 

h-1  - 1 » — - — — -  •  Z/-  (25b> 

{1  -  (l^X)  N2  (M2  +  (6M2  ^  )} 

r,  _  M  I1?'  ^  M  H)  -  ^  (H)1^2)-  (K)» 

{  1  -  (^)2  N2  (M2  +  (6M2  )  ) } 


(25c) 


21 


where 


2>  ■  K 


+  /sM2^>  +  /dM2 


Since 


M  =  M  +  <5M 
s 


We  have  by  squaring  and  averaging 


M2  =  M2  +  /SM2N 
s  \  / 


<6M2^  = 


2^  =  M2  -  M2 


When  the  material  is  isotropically  demagnetized 


(<}=  <aMz)=  7  <5mZ)  =  7  M' 


2 

V-'  -1  =  -  1  =  -4n  - ^ -  (32) 

.  N2  .  /  .2 

1  '  3  (UJM^ 

In  order  to  qualitatively  take  into  account  the  anisotropy  and  other 
internal  fields  that  were  dropped  from  Eq.  (22),  we  indicated  we  would  re¬ 
place  replace  6M  by  6M  +  <5H/N  in  our  final  results.  In  an  isotropically 
demagnetized  material  is  also  zero.  Thus  6H  is  parallel  to  <5M.  If 
we  define  Ha  as  ^|<5H|^,.Eq.  (32)  then  becomes 


(33) 


1=^  =  -IN 


(UM  +  y  Ha/N)7u 
1  -  ^  (uM  +  vHa/N)%2 


~T 


VI.  COMPARISON  WITH  EXPERIMENT 


The  next  question  to  be  considered  is  the  value  of  N.  If  our  do¬ 
mains  were  surrounded  by  nonmagnetic  material  the  sum  of  the  three  diag¬ 
onal  components  of  N  would  be  1  0  and  the  average  value,  N,  would  be  one 

3 

third.  However,  Polder  and  Smit  have  shown  that  for  some  domain  con¬ 
figurations  the  sum  of  the  diagonal  components  of  N  can  be  2.0.  Thus  we 
must  use  an  experimentally  determined  value.  In  Fig.  1,  we  have  p'  plotted 
vs  td.-Vu  for  several  materials  and  frequencies.  The  field  H  ,  a  com- 

1V1  cL 

bination  of  anisotropy  and  demagnetizing  fields,  is  not  well  known,  but  is 

much  less  than  4ttM  .  We  have  also  plotted  Eq.  (33)  for  H  =  0  and  4ttM  /10 

S  cl  s 

with  N  =  0.  60.  This  value  of  N  gives  the  best  fit  for  these  values  of  H  . 
With  larger  values  of  H  ,  a  smaller  value  of  N  would  be  required- 

3. 


We  thus  see  that  the  theory  gives  reasonable  agreement  with  the 
experimental  results  for  p'  vs  in  the  demagnetized  state.  We  shall 

now  examine  how  well  Eqs.  (25)  describe  the  dependence  of  p',  p1  ,  and  k' 
on  4ttM  as  one  cyles  about  the  B-H  hysteresis  loop.  First  one  should  note 
that  the  expressions  for  these  components  of  the  permeability  tensor  depend 


on 


(6M2), 


oM 


function  of 


x/ 
and  M 


and  <6M 


In  Eq.  (30)  we  showed  that  (6M  )  is  a  unique 


However,  its  components  are  not  unique  functions  of 


M  and  M  and  depend  in  addition  on  the  specific  domain  geometry  and  hence 
s 

on  hysteresis  history.  Thus  the  permeability  shows  hysteresis  as  a  function 
of  4ttM.  Two  states  having  the  same  magnetization  but  on  opposite  sides  of 
the  B-H  hysteresis  loop  have  different  permeability.  This  hysteresis  is 

4 

observed  experimentally.  Its  magnitude  and  even  sign  vary  from  material 
to  material.  Without  a  better  understanding  of  domain  dynamics  in  poly¬ 
crystalline  materials,  a  theoretical  description  of  how  (<5M2^>  and  ^6M2z^ 
vary  around  the  hysteresis  loop  is  not  possible. 
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in  the  Demagnetized  State 


We  next  note  the  difference  in  behavior  of  jx '  and  as  a  function 
of  4ttM.  The  prime  difference  is  the  term  M  H^)in  (jl1.  At  applied  dc  fields 


several  times  the  coercive  force  4ttM  is  still  well  below  saturation  but  is 
as  large  as  it  can  be  made  without  making  the  applied  dc  field  very  large. 

Since  does  not  go  to  zero  until  4ttM  =  4TrMg)  |P  and  \i'^  do  not  equal 

1. 0  at  this  point.  However,  (P  is  approaching  1 .  0  as  4itM  -*  4ttM  and 

/  Z  \  ^  s  4 

(  6M  >  -*0.  This  follows  from  Eq.  (Z5b)  and  is  observed  experimentally. 

On  the  other  hand  p'  is  observed  not  to  approach  1 . 0  as  4ttM  -»  4ttM  but 

s  4 

rather  flattens  out  and  eventually  turns  downward.  Earlier  explanations 
attributed  this  to  the  contribution  of  the  applied  dc  field  to  the  permeability. 
This  is  only  part  of  the  explanation.  The  actual  source  of  this  roll  off  in  jjl* 
vs  4ttM  is  the  term  M  .  While  the  applied  dc  field  contributes  to  <^H^)  , 
the  major  contribution  to  this  term  is  anisotropy  and  dc  demagnetizing  fields. 

It  should  be  noted  that  the  average  internal  effective  dc  field  taken  near 
saturation  and  H&  =■•  <^|  <5H  ,  the  average  magnitude  of  this  internal  effective 

dc  field  taken  in  the  demagnetized  state  are  related  quantities  but  not  neces¬ 
sarily  identical. 

To  -nderstand  in  detail  how  <^H)and  <^|<5H|^)  depend  on  anisotropy 
and  porosity,  it  is  necessary  to  have  a  better  understanding  of  domain  struc¬ 
ture  than  is  presently  available.  However,  it  can  be  seen  that  H  is  a  measure 
of  how  tightly  each  domain  is  bound  to  its  equilibrium  orientation.  Increasing 
porosity  puts  additional  constraints  on  the  domain  formation,  binding  the  do¬ 
mains  more  closely  to  their  equilibrium  position  and  increasing  ^H^  for 
partially  magnetized  states.  As  a  result  [P  will  not  come  as  close  to  1 .  0 
for  a  porous  materials  as  a  dense  one  with  the  same  magnetization  and  aniso¬ 
tropy.  It  is  observed  experimentally  that  microwave  phase  shifters,  such  as 
5 

PIP,  that  depend  on  the  change  in  |P  from  the  demagnetized  state  to  that  at 
B-H  loop  closure  for  their  phase  shift,  lose  phase  shift  rapidly  with  small 
amounts  of  porosity.  This  tends  to  confirm  the  previous  argument.  However, 
the  argument  can  be  rephrased  to  predict  the  H&  =  6H  also  increases 
with  porosity  and  hence  |(P  -  1 1  is  greater  for  porous  materials  than  dense 
materials  when  4ttM  =  0.  This  has  not  been  observed,  and  might  be  due  to 
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the  fact  that  the  demagnetized  state  is  more  random  and  hence  the  porosity 
may  have  less  effect  in  "stiffening"  the  domain  pattern  than  in  the  more  or¬ 
dered  state  near  B-H  loop  closure. 


The  remaining  prediction  of  this  theory  to  be  examined  is  the  de¬ 
pendence  of  k'  on  4rrM.  The  leading  term  of  Eq.  (26c)  is  the  same  as  Rado's 
result  as  is  expected.  The  corrections  include  terms  containing  <6M  \ 
and  (6M^)  and  hence  hysteresis  in  k'  vs  4ttM  is  predicted.  This  is  ob¬ 
served  experimentally.  However,  the  theory  has  this  hysteresis  going 

4 

to  zero  when  4ttM  =  0,  and  this  is  contradicted  by  experiment.  It  is  at 
this  point  that  the  deficiencies  of  the  model  become  apparent.  The  terms 
in  SH  which  we  neglected  as  uncorrelated  to  6M  do  have  a  measurable  ef¬ 
fect  although  small  compared  to  the  other  effects  predicted  by  the  theory. 

In  particular  by  ignoring  the  correlation  between  N  and  the  direction  of  6M 
we  lost  information  that  deferentiate j  the  two  states  with  4ttM  =  0  and  applied 
field  equal  to  ±  the  coercive  force. 


VII.  SUMMARY 


The  theory  presented  here  predicts  the  following  phenomena  which 
are  experimentally  observed: 


1.  [i1  and  fi'  differ  from  1  0  in  partially  magnetized  states, 

z 

2.  A  formula  is  given  for  the  dependence  of  |x'  and  |i'  in  the  de¬ 
magnetized  state  on  u.  /u  which  is  quantitatively  correct  for 

Vu  <  °- 9- 

3.  Hysteresis  in  p. ' ,  y.! ,  and  k'  vs  4ttM  is  to  be  expected.  The 
theory  allows  for  hysteresis  of  either  sign  and  does  not  predict 
its  magnitude, 

4.  approaches  1.  U  as  the  material  saturates, 

5.  (i1  does  not  approach  1 . 0  as  the  material  saturates  and  re¬ 
mains  below  it  by  an  amount  given  by  an  effective  internal  dc 
field  made  up  of  anisotropy,  demagnetizing,  and  the  applied 
field. 


The  theory  also  predicts  that  hysteresis  in  k'  goes  to  zero  as 
4irM  goes  to  zero.  This  is  contradicted  by  experiment. 
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This  appendix  contains  the  algebra  necessary  to  obtain  the  perme¬ 
ability  tensor  components  as  given  in  Eq.  (25)  from  the  formula  for  6K  given 
in  Eq.  (27). 

As  a  notational  convention  A  X  B  XC  will  mean  A  X  (B  XC).  We 
make  use  of  the  following  vector  identities: 


A  X  B  XC  =  B(A  •  C)  -  C(A  •  B)  (A  1 ) 

A  •  B  X  B  =  0  (A2) 

and  the  definitions 

M  =  M  +  6M  (A3) 

s 

<6M2)  =  +  <5M2)+  (<5M2)  (A4) 

and  the  symmetry  condition 


(6M2)=  (6M2)  .  (A5) 

We  begin  by  twice  taking  the  cross  product  from  the  left  of  Eq.  (24) 

with  M  • 
s 

M  XSR  =  {^M){M  X  6lx(h)  +  I  XI  X  6H  -  M  X  (sm  X  6h) }  (A6) 

M  X  M  X  <5h  =  (in?)  {M  X  M  X  <5M  x/h\  +  M  X  M  X  6h  -  M  XI 
ss  ju  '  s  s  \/ss  ss 

x(<5MX<5h)}  .  (A7) 
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Making  use  of  Eqs .  (Al),  (A2),  and  (A3),  Eq.  (A7)  becomes 

M  XM  X  =  M  •  (M  x(fi\  )  -  M2  (6Mx(fi\)  -  M2  (M  X  6fi) 

-  MgMg  •  (m  X  6h)  +  M2  (fiM  X  afi)  }  .  (A8) 

Substituting  Eq.  (A8)  into  Eq.  (A6) 

M  X  6H  =  [M  X  6Mx(h\+  {-MM  •  (MX  (fi\)  -  M2  (6MX  (fi  )) 

S  JCJ  S  \  /  JGJ  S  S  \  /  S  '  ' 

-  M2  (MgX  6E)  -  MgMs  •  ^5MX6fi^+M2  <6MX  6h^}  -  MgX  <^M  X  6H^)  ] 

(A9) 

2 

Taking  the  spatial  average  of  both  sides,  noting  that  M  is  a  constant  gives  « 

s 

<Msx  6K)  =  [<ra3XS“x<K»+  T('^  ' (Mx  <h>> 

-M^  /fjMz(K))  (msX  6H)-{MsHs)  •  4mx  oH) 

+  M2  (fiM  X  5h)  }  -  (Ms)  X  (m  Xafi)  }  (A  10). 

Using  Eq.  (A3)  and  noting  that  the  product  of  constants  and  one  spa¬ 
tially  fluctuating  variable  has  an  average  of  zero,  we  get,  after  collecting 
terms  in  Eq.  (A10) 

(fiM  X  6h)  =  [  (m  X  5M  ><(  h))  +  {-<MgMg)-  (M  X  (in) ) 

-  (msms)  •  (m  x  ah^ }  -  m  x  (m  x  an)  ]  .  (Ai  i) 

By  symmetry 


<3M  6M  \  =  <6M  3M  \  =  (6M  6M  \  =  0 
\  x  y/  \  y  z/  \  z  x/ 


(A  12) 


Therefore 


/M  M  V  ~K  =  /M2  \  A  x  +  /M2  \ 
\  s  s  /  \  sx  /  x  \  sy/ 

where  A  is  any  constant  vector, 


/\ 

A  y  + 
y 


<mL)  ■  <K> 

<Msy>  =  <SMy> 

(Msz)  =  mZ  +  (6Mz) 

^  /\  A 

and  x,  y,  and  z  are  unit  vectors. 


(A13) 


(A  14) 


In  addition 

<6M  X  «5M  x(E))=  (<5M<5m)  •  /h)-  <6M2)  (h) 

=  «K)  -  (5m2>  <hx>  +  (<6My>  -  <6m2>  (hy)  y 
+  (  (fiM2)  -  (6M2))  (hz)z  .  (A 1 5) 

If  we  write  Eq.  (All)  in  component  form,  making  use  of  Eqs.  (A13)  and  (A  15), 
we  have  three  simultaneous  equations  for  the  three  components  of  (6M  X6h^  . 
Without  loss  of  generality  we  can  again  take^h^)  to  be  in  the  x-z  plane  so  that 

(hy>=  0 

<«  X  «K)X  =  [  (  <«M^)  -  («M2))  <V>+  <«M  X  ahp 

+  M  (m  X  6E)y] 

(dm  X  6h)  =  {-(6M2)m  (h^  -  (6M2^)  (dM  X  6h)  }  -  M  (M  X  6E)x] 

v  3  ^  3  ^  y  y  y 

(aMXS%  -  ±^Il(«Mz2>-(«M2))<hz)+i^{-(M2+  (6M2))  <6BX6E)z}] 

(A  16) 
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Collecting  terms  and  making  use  of  Eq.  (A5)  gives 


(6M2)}  (SmX 


)  M  <610  X  6H)x  +{1  +  }(6H  X  ST5>y  H±$?>  («M2)  M  (hjt> 


{  1  ♦  (il^i )  (M2  t  (6Mz2))}<6M>  X  SR\  =  (  (6M2)-  (SM2))  (hj 


(A  17) 


Solving  these  equations  yields 


i 


i 

i 


(A  18) 


Since  we  are  concerned  only  with  the  region  u^/u;  <  1 ,  all  terms  of 
the  form  (47ryN/uj)^  (^6M^  or  or  <^>M^  )  will  be  much  less  than  1.  0. 

We  thus  drop  all  but  the  most  significant  terms  in  (A  18)  and  make  use  of 
Eqs.  (A4)  and  (A5)  leaving 


<6M  X  <5fi)x  = 


(A  19a) 
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MICROWAVE  BEHAVIOR  OF  PARTIALLY  MAGNETIZED  FERRITES 

by 

Ernst  Schlomann 
Raytheon  Research  Division 
Waltham,  Massachusetts  02154 


ABSTRACT 


The  high  frequency  permeability  of  partially  magnetized  ferrites  is 
calculated  for  some  simple  domain  configurations,  comprising  only  "upl¬ 
and  ,,down"-domains.  The  method  used  is  based  upon  the  magnetostatic 
approximation  and  neglects  exchange  effects,  but  is  otherwise  substan¬ 
tially  rigorous.  The  components  of  the  effective  permeability  tensor 
(ratio  of  average  induction  to  average  magnetic  field)  in  general  depends 
upon  details  of  the  domain  configuration  in  addition  to  the  net  dc  magnetiza¬ 
tion.  When  the  dc  magnetization  is  cycled  between  the  two  states  of  complete 
magnetization  the  high  frequency  permeability,  considered  as  a  function  of 
the  dc  magnetization,  in  general  shows  hysteresis.  Detailed  calculations 
of  the  high  frequency  permeability  have  been  carried  out  for  the  case  in 
which  the  domain  configuration  is  cylindrically  symmetric,  i.  e.  ,  invariant 
under  rotation  around  the  direction  of  magnetization.  For  any  such  domain 
configuration  the  two  relevant  components  and  of  the  effective  per¬ 
meability  tensor  obey  the  relation  p2^  -  k2^  =  const,  regardless  of  the 
number  of  domains  and  their  relative  size.  This  general  theorem  allows  a 
simple  derivation  of  the  (isotropic)  permeability  in  the  completely  demag¬ 
netized  state,  giving 


A 


eff 


=  f  {  [  (w/y)2  -  (Ha  +  4uM0)2]/l  (W/Y)2  -  H2]H 


where  u  is  the  frequency,  H  the  anisotropy  field  and  Mn  the  saturation 

3.  U 

magnetization. 

*This  research  was  supported  in  part  by  the  Advanced  Research  Projects 
Agency  of  the  Department  of  Defense  and  was  monitored  by  Rome  Air 
Development  Center  under  Contract  No.  F30602-69-C-0026. 
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1.  INTRODUCTION 


In  many  of  the  applications  of  ferrites  in  microwave  devices  such 
as  phaseshifters,  switches  and  rotators  the  magnetic  material  is  only 
partially  magnetized.  This  mode  of  operation  has  the  advantage  that  the 
microwave  properties  of  the  material  can  be  changed  significantly  by 
the  application  of  a  rather  modest  magnetic  field. 

Because  of  these  practical  applications  the  microwave  be  ha',  lor  of 

partially  magnetized  ferrites  has  been  the  subject  of  considerable  research 

1  2 

from  the  early  days  of  the  ferrite  microwave  technology  ’  uni-]  very 
3  “  6 

recently.  It  has  been  found  that  the  off-  diagonal  components  of  the  sus  ¬ 
ceptibility  tensor  are  substantially  proportional  to  the  net  magnetization, 
and  that  the  diagonal  components  are  negative  ard  increase  in  magnitude 
with  increasing  saturation  magnetization  and  decreasing  frequency.  It  has 
also  been  known  for  a  long  time  that  the  partially  magnetized  state  is 
characterized  by  very  high  losses,  when  the  frequency  ic-  is  less  than,  or 
of  the  order  of,  =  y4ttM0  (y  =  gyro  magnetic  ratio,  MQ  =  saturation 
magnetization). 

The  first  successful  theory  of  the  microwave  behavior  of  partially 

7  8 

magnetized  ferrites  is  due  to  Rado.  ’  By  suitably  averaging  the  magnetic. 

equation  of  motion  and  using  certain  simplifying  assumptions  concerning 

the  correlation  of  magnetization  and  magnetic  field,  he  deduced  that  the 

off-diagonal  component  of  the  susceptibility  tensor  should  be  giver,  by 

YM/io  where  M  is  the  net  magnetization.  This  result  is  in  quite  good 

agreement  with  the  experimental  observations.  According  to  RadoEs  theory 

the  diagonal  components  of  the  susceptibility  tensor  would  be  expected  to 

vanish  (for  zero  dc  magnetic  field  ir.  the  limit  of  zero  loss';.  Similar 

9  10 

results  have  also  been  derived  by  Soohoo  and  by  Allen.  A  more  refired 
theory,  based  upon  a  similar  approach  but  taking  the  correlation  of  local, 
magnetic  field  and  local  magnetization  more  carefully  into  account  has 
recently  been  proposed  by  Sandy,  *  This  theory  yields  finite  diagonal 
components  and  agrees  quite  well  with  the  experimental  data,  provided 
that  a  certain  parameter  which  occurs  in  the  theory  (an  effective  rf  demag¬ 
netizing  factor)  is  properly  chosen. 
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A  very  important  suggestion  concerning  the  microwave  behavior 

12 

of  partially  magnetized  ferrites  has  been  made  by  Polder  and  Smit. 

They  pointed  out  that  in  a  sample  containing  oppositely  magnetized 
domains  the  resonant  frequencies  of  the  normal  modes  of  vibration 
extend  upwards  to  the  vicinity  of  ^ttMq.  The  "low  field  loss"  generally 
observed  in  ferrites  when  the  frequency  is  sufficiently  low  can  be 
attributed  to  these  Polder -Smit  resonances. 

Various  attempts  have  been  made  to  develop  a  quantitative  theory 
of  the  microwave  behavior  of  partially  magnetized  ferrites  on  the  basis  of 
the  Polder -Smit  suggestion.^’  *  These  calculations  are  based  upon 
some  ad  hoc  assumptions  (for  instance,  concerning  the  distribution  of 
resonant  frequencies)  and  are,  therefore,  of  limited  use,  even  though 
good  agreement  with  the  experimental  data  has  been  achieved  in  some  cases. 

In  the  present  paper  a  new  theory  of  the  microwave  behavior  of 
partially  magnetized  ferrites  is  developed.  This  theory  is  in  a  certain 
sense  complementary  to  those  previously  proposed,  since  it  avoids  many 
of  the  shortcomings  of  the  previous  theories  but  also  has  some  shortcomings 
of  its  own.  The  microwave  behavior  (at  low  power  levels)  can  generally 
be  characterized  by  an  "effective"  susceptibility  tensor,  v/hich  relates  the 
spatial  average  of  the  rf  magnetization  to  the  spatial  average  of  the  rf 
magnetic  field.  The  aim  of  the  theory  described  in  the  present  paper  is 
to  calculate  this  effective  susceptibility.  From  the  magnetic  equations  of 
motion  one  can  easily  determine  the  "local"  or  "microscopic"  susceptibility 
at  any  point  inside  the  sample,  if  the  direction  of  the  local  saturation 
magnetization  at  this  point  is  known.  The  microscopic  susceptibility 
relates  the  local  rf  magnetization  to  the  local  rf  magnetic  field.  In  order 
to  calculate  the  effective  susceptibility  it  is  necessary  to  determine  the 
local  rf  magnetic  field  first.  In  some  of  the  previous  theories  this  problem 
is  circumvented  by  postulating  a  relationship  between  the  local  magnetization 
and  the  local  magnetic  field  (an  effective  rf  demagnetizing  factor).  Such  a 
procedure  can  not  be  rigorous,  because  the  rf  magnetic  field  at  any  point 


inside  the  material  is  actually  not  determined  by  the  rf  magnetization  at 
that  point,  but  rather  by  the  distribution  of  rf  magnetization  throughout 
the  sample.  In  the  present  theory  the  calculation  of  the  local  rf  magnetic 
field  is  considered  to  be  the  central  issue.  Once  this  problem  is  solved 
the  effective  susceptibility  can  easily  be  obtained. 

Because  of  the  complexity  of  the  problem  at  hand  the  calculation 
described  above  can  only  be  carried  out  for  some  highly  idealized  domain 
configurations.  In  the  present  paper  we  confine  attention  to  the  case  in 
which  the  sample  contains  only  two  types  of  domains,  which  are  aligned 
either  parallel  ("up -domains'*)  or  antiparallel  ("down -domains")  to  some 
reference  axis  (which  is  taken  as  the  z-axis  of  the  coordinate  system).  It 
is  shown  in  Section  2  that  for  this  case  the  problem  of  finding  the  local  rf 
magnetic  field  reduces  to  solving  a  simple  differential  equation,  subject  to 
certain  (not  so  simple)  boundary  conditions  at  the  domain  boundaries.  A 
complete  solution  of  the  problem  is  constructed  in  Section  3  for  the  case 
in  which  the  sample  is  a  circular  cylinder  and  contains  a  cylindrically 
symmetric  domain  structure.  It  is  shown  that  for  any  domain  configuration 
of  this  type  the  components  of  the  effective  permeability  are  related  by  a 
general  theorem  (Section  3a).  Explicit  results  concerning  the  variation  of 
the  effective  susceptibility  with  the  net  dc  magnetization  are  derived  for 
two  simple  examples;  the  case  of  three  domains  of  arbitrary  volume 
(Section  3b),  and  the  case  in  which  all  comains  of  the  same  type  have  equal 
volume  (Section  3c).  The  results  are  discussed  and  critically  reviewed  in 
Section  4. 


2.  GENERAL  FORMULATION  OF  THE  THEORY 


As  previously  mentioned  we  assume  that  all  domains  are  oriented 
either  parallel  or  antiparallel  to  the  z-axis.  This  represents  a  stable 
domain  configuration  provided  that  the  domains  are  cylinders  whose  axes 
are  aligned  with  the  z-axis.  The  cylinders  can  have  arbitrary  cross- 
sections  in  the  x-y  plane  as  shown  in  Fig.  1.  The  domains  in  which  the 
dc  magnetization  is  parallel  to  the  positive  z-axis  (up-domains)  are  denoted 
by  the  subscript  u,  those  in  which  it  is  antiparallel  to  the  positive  z-axis 
(down-domains)  by  the  subscript  d. 

The  local  permeability  is  the  same  throughout  each  type  of  domain, 
but  differs  from  one  type  of  domain  to  the  other.  We  characterize  the  x 
and  y  components  of  the  local  permeability  tensors  for  the  two  types  of 
domains  by 


U  “  IK 

ru  u 


IK  U 

U  rU 


(1) 


The  tensor  components  p  and  k  can  easily  be  calculated  from  the  magnetic 
equations  of  motion.  If  loss  is  neglected,  one  finds 


u,  d 


=  1  + 


H  .4TrMn 
u,  d  0 


4ttMq  u/y 

Cu,d  ~  u2  .  /  .2 

Hu,d"  (u/V> 


(2) 


where  Hu  ^  is  the  effective  dc  magnetic  fieldstrength  in  the  direction  of  the 
(local)  dc  magnetization.  If  the  material  has  uniaxial  crystalline  anisotropy 
the  effective  dc  magnetic  field  arising  from  this  source  would  be  the  same 


in  both  types  of  domains.  An  externally  applied  dc  magnetic  field,  if  it  is 
not  compensated  by  demagnetizing  fields,  would  add  to  the  effective  dc  field 
in  one  type  of  domain  and  subtract  from  that  in  the  other.  Actually  the 
applied  field  is  almost  entirely  compensated  by  demagnetizing  fields.  By 
way  of  approximation  it  is,  therefore,  permissible  to  set  pu  =  p^,  =  k  . 

In  the  following  we  adopt  this  approximation  for  the  sake  of  simplicity. 
Furthermore,  the  anisotropy  field  is  usually  much  smaller  than  to/y.  Under 
these  conditions  Eq.  (2)  can  be  approximated  by 

P  =  1  .  x  =  uM/u  (3) 

where  uM  =  y4ttMq. 

The  effect  of  dissipation  upon  the  local  permeability  can  be  taken 
into  account  by  adding  a  phenomenon  logical  damping  term  to  the  equations 
of  motion,  p  and  k  are  then  in  general  complex.  In  the  limit  of  zero 
effective  magnetic  field  the  tensor  components  are  now 

P  a  1  -  ifi  .  x  =  uM/u  .  (4) 

The  off-diagonal  tensor  component  k  is  real  in  this  limit.  6  is  related 

1  4 

to  the  "effective"  linewidth  AHg^  that  characterizes  the  losses  of  the 
saturated  material  in  the  region  of  very  small  dc  fields  by 

6  -  Y2  2ttMq  AHeff/co2  .  (5) 

If  the  profile  of  the  resonance  line  is  Lorentzian  the  effective  linewidth  that 
occurs  in  Eqs.  (4)  and  (5)  equals  the  width  AH  of  the  resonance  peak.  In 
many  polycrystals,  however,  AHg^  is  much  smaller  than  AH. 

Consider  now  the  rf  magnetic  field  h(:r )  associated  with  some  arbitrary 
distribution  of  the  rf  magnetization  QjOJ.  We  calculate  this  field  subject  to 
the  quasistatic  approximation.  Thus  jj(£}  can  be  represented  as  the  gradient 
of  a  potential  \&(:r  ) 


Since  the  divergence  of  b  =  (jl  •  h  vanishes  the  potential  ^  satisfies  the 
differential  equation 


which  is  well-known  from  the  theory  of  magnetostatic  modes.  ^  It 
should  be  noticed  that  the  same  equation  (7)  applies  in  both  types  of  domains. 


In  the  present  context  we  need  to  consider  only  such  field  distribu¬ 
tions  that  are  independent  of  z.  In  this  case  the  differential  equation  that 
determines  itr  reduces  to  the  two-dimensional  Laplace  equation 


The  boundary  conditions  to  which  the  potential  is  subject  arise  from 
the  requirement  that  the  tangential  component  of  the  magnetic  field  h^  and 
the  normal  component  of  the  magnetic  induction  must  be  continuous.  Let 
3^/3  s  denote  the  tangential  derivative  and  3^/3  n  the  normal  derivative  of 
the  potential  taken  at  the  domain  boundary.  The  potential  ^  and  ^  (applicable 
inside  the  up-  and  down-domains  respectively)  then  has  to  satisfy  the  condi¬ 
tions 


3^  3tf, 
_ u  _  _ d 

3  s  3  s 


/3¥ 

u 

3* 

u 

3  n 

i 

3  n 

=  ix 

i8s 

, 

+  3s 

i 

(9) 


Here  it  is  assumed  that  n,  s,  z  constitutes  a  right-handed  coordinate  system 
(as  does  x,  y,  z). 


T 


In  the  following  it  will  be  assumed  that  the  sample  is  surrounded 
by  air.  In  the  exterior  of  the  sample  the  potential  SH ,  therefore,  also 
satisfies  the  Laplace  equation  (8).  At  the  sample  boundary  the  tangential 
component  of  h^  and  the  normal  components  of  ^  are  continuous. 

Assume  now  that  the  differential  equation  for  the  potential  has  been 
solved  subject  to  all  the  boundary  conditions  at  the  domain  boundaries  and 
the  sample  surface.  How  is  this  solution  related  to  the  microwave  behavior 
that  is  directly  observable? 

The  magnetic  response  of  a  small  sample  of  ellipsoidal  shape  can 
be  characterized  by  the  "internal"  or  the  "external"  effective  susceptibility. 
The  former  relates  the  average  rf  magnetization  to  the  average  rf  magnetic 
field  inside  the  sample,  the  latter  relates  it  to  the  rf  magnetic  field  existing 
outside  the  sample.  For  experiments  with  small  samples  the  external  sus¬ 
ceptibility  is  more  closely  related  to  the  directly  observed  physical  quantities 
(e.  g. ,  change  of  frequency  and  l/Q  of  a  resonant  cavity).  The  ir.terna.1. 
susceptibility,  on  the  other  hand,  is  a  more  fundamental  quantity,  because 
it  is  not  dependent  upon  the  shape  of  the  sample.  We  shall,  therefore,  state 
our  results  in  terms  of  the  internal  effective  susceptibility.  It  is  well  known 
that  the  two  types  of  susceptibilities  are  related  by 

y  =  (1  +  4tt  y.  N.)  ^  y.  (10) 

where  Nj.  is  the  transverse  demagnetizing  factor.  All  susceptibilities  and 
permeabilities  mentioned  in  the  remainder  of  the  paper  are  "internal".  The 
subscript  "in"  will  henceforth  be  omitted. 

The  microwave  behavior  of  partially  magnetized  ferrites  can  most 
economically  be  characterized  by  the  "circular"  susceptibility,  i.  e.  ,  the 
response  to  a  rotating  rf  magnetic,  field.  If  the  domain  structure  does  not 
distinguish  any  direction  in  the  x-y  plane  on  a  macroscopic  scale  the  sus¬ 
ceptibility  tensor  must  be  diagonal  with  respect  to  such  a  rotating  magnetic 
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field  We  calculate  primarily  the  effective  circular  permeability 
p.  e^{u)  f°r  Positive  sense  of  rotation.  The  corresponding  permea¬ 
bility  for  the  negative  sense  of  rotation  is  then  obtained  by  inverting  the 
sign  of  the  frequency  and  taking  the  complex  conjugate 


(«)  =  *  (-u) 

eff  eff 


The  effective  permeability  tensor  has  the  same  structure  as  the  local 
permeability  tensor  [given  by  Eq.  (1)]  except  that  the  components  pe^  and 
are  now  defined  by 


eff  =  ^eff  T  Keff 


In  the  next  section  the  calculation  of  the  effective  permeability  is 
explicitly  carried  out  for  a  simple  domain  structure. 


3.  APPLICATION  TO  SAMPLES  WITH  CYLINDRICAL  DOMAIN  STRUCTURE 


3a.  General  Discussion 


Consider  now  a  sample  in  the  shape  of  a  circular  cylinder,  with  a 
cylindrically  symmetric  domain  structure  such  as  shown  in  Fig.  2.  We 
assume  (arbitrarily)  that  the  domain  at  the  core  of  the  cylinder  is  a  down- 
domain.  The  significance  of  this  assumption  is  further  discussed  later  in 
this  section.  The  outer  radii  of  the  domains  are  denoted  by  r^(n  =1,  2. .  .N). 

The  fundamental  solutions  of  Laplace1  s  equation  in  cylindrical  coor¬ 
dinates  are 

T  m  -im<t>  -m  -imd)  ,n, 

=  r  e  T  ,  r  e  Y  (13) 


where  m  is  a  positive  integer.  Consider  the  first  of  these  solutions  for 
m  s  1,  j.  e. ,  s  r  e  =  x  -  iy.  Since  the  time  dependence  is  assumed  to 
be  exp  (Lot)  it  follows  that  the  magnetic  field  associated  with  this  potential  is 


,  „  not  . 

h  =  Re  e  =  cos  tot 
x 


h  =  Re  (  i 

y 


sin  tot 


(14) 


Thus  for  m  =  1  the  first  potential  :n  Eq.  (13)  describes  a  rotating,  uniform 
magnetic  field.  The  sense  of  rotation  depends  upon  the  sign  of  to,  being 
positive  for  positive  to.  In  addition  to  the  solutions  given  by  Eq.  (13)  the 
Laplace  equation  also  admits  solution  of  the  form  rm  elm^  and  r  m  e'*'111^. 

The  fields  associated  with  these  pole-' vials  are  identical  with  these  of  the 
potentials  given  by  Eq.  (j3)  provided  that  the  sign  of  to  is  reversed.  Thus 
these  additional  solutions  need  not  be  considered  separately  in  the  calculation 
if  to  is  allowed  negative  as  well  as  positive  values. 
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In  general  the  potential  inside  each  domain  as  well  as  in  the  exterior 
of  the  sample  can  be  expanded  in  terms  of  the  fundamental  solutions  given 
by  Eq.  (13).  In  the  present  context  we  can  confine  attention  to  those  poten¬ 
tials  that  describe  an  rf  magnetic  field  which  reduces  to  a  homogeneous 
rotating  field  at  points  far  away  from  the  sample.  For  r-*  °o  the  potential 
must  be  of  the  form 

^  (r  -*  °o)  =  hQr  e  ^  (]  5) 

where  hQ  is  the  amplitude  of  the  rotating  magnetic  field.  Because  of  the 

orthogonality  of  the  exponential  functions  for  different  m,  the  expansion  of 

the  potential  inside  the  sample  then  also  contains  only  the  terms  corresponding 

th 

to  m  =  1.  Thus  the  potential  for  the  n  domain  is 

Jj  =  (A  r  +  B  r_1)  e'1^  (16) 

where  An  and  Bn  are  arbitrary  coefficients,  to  be  determined  by  the  boundary 
condition.  For  the  core -domain  (n  =  1) 


B 


(H) 


to  avoid  a  singularity  at  r  *  0.  For  the  region  outside  the  sample  (n  =  N  +  1) 


AN+1 


=  h 


0 


to  conform  with  Eq.  (15). 


(18) 


The  radial  and  circumferential  components  of  h  and  b  inside  the 

th  ^  ^ 

n.  domain  are  according  to  Eqs.  (16),  (6)  and  (1) 
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P  f  j Byr  ■  tl  jli 


Equations  (20-22)  constitute  an  inhomogeneous  set  of  2N  linear 
equations  for  the  2N  arbitrary  coefficients  An  and  B^.  Each  of  these  coef¬ 
ficients  is,  therefore,  uniquely  determined.  Not  surprisingly,  each  is 
proportional  to  hQ. 
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Consider  now  the  spatial  average  of  the  rf  magnetic  field  and  the 

rf  magnetic  induction.  The  x-  and  y-  components  of  the  magnetic  field 
th 

inside  the  n  domain  are  according  to  Fig.  2  and  Eq.  (19) 

h  =  h  cos  <{>  -  h ,  sin  <b  =  A  -  B  r^e  ^ 
xn  rn  T  <}>n  T  n  n 

h  =  h  sin  <j)  +  h ,  cos  <j)  =  -i  [A  +  B  r  2  e  2^]  .  (23) 

The  contribution  of  the  B -coefficients  vanishes  when  averaged  o\  er  <{>. 
It  can  be  shown  by  substantially  the  same  reasoning  that  this  is  also  true  for 
the  rf  magnetic  induction.  Thus  only  the  A -coefficients  contribute  to  the 
average  magnetic  field  and  to  the  average  magnetic  induction. 


In  order  to  obtain  the  effective  circular  permeability  we  have  to  take 
the  ratio  of  the  spacial  averages  of  the  magnetic  induction  and  the  magnetic 
field.  Thus 


H+  eff 


(p+K)  r2  Aj  +  (p-K)  (r2  -  r 2)  A2  +  (p+«)  (r2  -  r2)  A 3  + 
r(  A!  +  [t\  -  r()  A2  +  <r*  -  t\)  Aj  +  . . . 


(24) 


For  further  progress  the  coefficients  A^  must  be  determined.  For 
convenience  we  define  a  ’’vector” 


* 


n 


\  B 


n 


Equation  (21)  is  then  equivalent  to 


A  ,  ,  =  M  A 
~n+l  ~n  ~n 


where  the  matrix  is  given  by 


(25) 


(26) 
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I 


M 


-n 


/  1  T  tr  T  or  \ 
'  n  ' 


>  ±  cr  r  1  ±  <r 
n 


(27) 


(upper  sign  for  even  n,  lower  for  odd  n  as  before).  Here 


=  k/M- 


(28) 


Equation  (26)  holds  also  for  n  =  1  except  that  in  this  case  Bj  =  0  [compare 
Eq.  (20)].  Using  Eq.  (26)  all  coefficients  An(n  =  2,  3. ..)  can  readily  be 
expressed  in  terms  of  Aj. 

It  is  apparent  from  the  preceding  results  [in  particular  Eq.  (24)] 
that  the  effective  permeability  in  general  depends  upon  details  of  the  domain 
configuration  (such  as  the  number  of  domains  and  their  size)  as  well  as  upon 
gross  features  (such  as  the  net  dc  magnetization).  For  any  reasonably 
complex  domain  structure  the  theory,  therefore,  yields  a  result  which 
depends  upon  a  vast  number  of  adjustable  parameters.  This  may  facilitate 
obtaining  a  good  fit  to  experimental  data,  but  the  significance  of  such  a  fit 
is  doubtful  unless  it  can  be  shown  that  the  domain  structure  actually  present 
is  similar  to  that  assumed  in  the  calculation. 


It  is  fortunate,  therefore,  that  certain  results  of  the  theory  are  not 
dependent  upon  the  details  of  the  domain  configuration.  It  is  shown  in 
Appendix  1  that  for  any  cylindrically  symmetric,  but  otherwise  arbitrary 
domain  configuration 

^eff  "  Keff  =  11  "  K  ‘  ^ 

This  theorem  has  several  intei  esting  consequences.  It  is  known  from  Rado’ s 
7  8 

work  ’  that  is  substantially  proportional  to  the  net  magnetization.  (The 
present  theory  gives  similar  results,  as  further  discussed  in  Sections  3b  and 
3c..)  Since  the  right  hand  side  of  Eq.  (29)  is  independent  of  the  domain 


configuration,  it  follows  that  considered  as  a  function  of  the  net 
magnetization  goes  through  a  minimum  at  that  magnetization  at  which 
Kg££  vanishes.  The  minimum  value  of  is  independent  of  the  domain 
configuration  being  given  [according  to  Eqs.  (29)  and  (2)]  by 


^eff 


min 


(10/ y) 2  -  (Ha  +  4ttM0)2 


(w./y)' 


H 


{30) 


Here  it  has  been  assumed  that  the  effective  magnetic  field  inside  the 


domains  is  the  anisotropy  field  H  .  For  a  sample  containing  many  very 

cl 


small  domains  the  state  of  vanishing  is  very  close  to  that  of  vanishing 
net  magnetization.  Thus  the  effective  permeability  given  by  Eq.  (30)  can 
be  identified  with  that  applicable  in  the  completely  demagnetized  state.  It 


should  be  noticed  that  it  is  real  for  u/y  >  +  4irM0  but  imaginary  for 


H  <  —  <  H  +  4trM,..  With  suitable  corrections  to  allow  for  the  random 

a  y  a  0 

orientation,  of  the  domains  and  for  the  finite  linewidth  of  the  material  Eq.  (30) 


has  been  found  to  agree  quite  well  with  the  observed  behavior  of  the  real  and 
imaginary  parts  of  in  the  demagnetized  state. 


3b.  Application  to  Samples  Containing  Only  Three  Domains 


In  this  section  and  the  one  following  ?t  the  general  results  derived 
in  the  preceding  section  are  explicitly  evaluated  for  two  special  cases.  The 
results  given  in  Section  3b  apply  to  a  sample  containing  only  three  domains 
(of  arbitrary  volume),  those  given  in  Section  3c  to  a  sample  containing  an 
even  number  N  (N  =  2,  4,6,  .10,  20,  40,  80)  of  domains  provided  that  domains 
of  the  same  type  have  equal  volume.  In  both  of  these  cases  the  results 
depend  upon  a  single  "adjustable"  parameter  as  well  as  upon  the  net.  magnet¬ 
ization.  In  the  three -domain  case  the  adjustable  parameter  is  the  volume 
ratio  of  the  two  down-domain;  in  the  many-domain  case  this  parameter  is 
the  total  number  (N)  of  domains,  According  to  Eqs.  (25)  and  (2b) 


A2  =  (1  +  tr)  A1 


A3  =  [  1  -  cr2  +  cr2  (r1/r2)2]  Aj 


The  relative  dc  magnetization  (i.  e.  ,  the  volume  difference  of  up 
domains  and  down-domains  divided  by  the  total  volume)  is 


,  „  2  ,2  1,1 
i  -  (2r2  -  2rj  -  rj/r^ 


and  the  volume  ratio  of  the  two  down-domains 


.2  l.,l 
p  =  (r3  -  r  )/r 


After  trivial  calculations  one  obtains  from  Eqs.  (24)  and  (31-33) 

,  1H  2  P(i  -  f) 


2[2  +  p(l  +  &)] 

Consider  now  the  real  parts  of  the  components  of  the  internal 
effective  susceptibility  tensor.  As  long  as  <  0.  9  the  losses  have 

only  a  very  small  influence  on  these  components,  and  are  therefore 
neglected  in  this  discussion.  According  to  Eqs.  (34),  (11),  (12)  and  (3) 
taking  tr  as  real 


=  1  -K  -  &  t 1  -  ^  p  zTp  T&T 1  /w 


e«  -r(, ■  -  Ml-i  )  (1-4)  ,1+„2,  _  ,  rp(l-n  . 

P  U  4  1  2  +  p  (1H)  '  y  J  2[ 2  +  p(lH)  ]  J 


(36) 


If  Ha  «  u/y  we  may  approximate  tr  by  w^/u).  It  may  be  seen  that  for 
tr  «  1  the  leading  term  in  k1^  is  cr£  =  y4ttM/u)  in  agreement  with  Rado's 
result.  It  should  be  noticed,  however,  that  K*e^  as  given  by  Eq.  (35)  does 
not  necessarily  vanish  for  t,  =  0. 


In  Figs.  3  and  4  the  tensor  components  p'e^  and  k*  ^  are  plotted 

versus  relative  magnetization  L  for  w^/u)  =  0.  9  and  different  fixed  values 

of  the  parameter  p  assuming  H  =0.  For  comparison  the  result  of  Rado's 

theory  is  also  shown  in  Fig.  4  as  a  broken  line.  It  may  be  seen  that  for 

p  =  4  the  versus  t,  curve  is  approximately  symmetric,  the  k'  .  versus 

t,  curve  approximately  antisymmetric  in  £,.  In  this  connection  it  is  important 

to  remember  that  the  calculation  as  presented  is  based  upon  the  assumption 

that  the  core -domain  is  magnetized  "down".  The  curves  applicable  for  the 

case  in  which  the  core -domain  is  magnetized  "up"  can  be  obtained  by  reflecting 

the  p-curves  around  the  vertical  axis  and  the  k  -curves  around  origin.  When 

the  static  magnetization  is  cycled  between  the  completely  magnetized  states 

the  core-domain  will  most  likely  reverse  its  direction  at  each  half  cycle. 

The  effective  permeabilities  will,  therefore,  exhibit  a  hysteresis  when 

plotted  as  a  function  of  the  static  magnetization.  Such  hysteresis  effects 

5 

are  indeed  observed  experimentally. 


.  The  effect  of  a  finite  anisotropy  field  upon  the  curves  shown  in  Figs.  3 
and  4  can  easily  be  taken  into  account.  In  Eqs.  (35)  and  (36)  the  anisotropy 
field  enters  only  through  cr  and  p.  Since  the  curves  are  valid  for  a  given 


c r,  they  apply  for  any  combination  of  4trMn,  u,  and  H  for  which 
2  2  u  a 

c r  =  k/p  =  Wjyjw/[  u  -  y  ^a(Ha  +  4ttMq)]  has  the  appropriate  value.  The 

ordinate  then  represents  Peff/p  an^  k eff/ ^  rather  than  pg^  and  itself. 

The  curves  shown  in  Figs.  3  and  4  apply  for  instance  also  when  =  0.  817, 


YHa/uJ  =  0.  1. 


The  value  of  pg^  in  the  completely  magnetized  states  (i.  e.  , 


i  =  ±  1)  is  smaller  than  unity  when  the  anisotropy  field  is  finite. 


13:1 
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Consider  now  the  imaginary  parts  of  the  effective  circular 
permeabilities.  For  simplicity  we  confine  attention  to  the  case  in  which 
the  anisotropy  field  is  negligible,  i.  e.  ,  the  microscopic  permeability  is 
given  by  Eq.  (4).  According  to  Eqs.  (34),  (11),  and  (4)  we  can  generally 
express  the  imaginary  parts  of  the  effective  circular  permeabilities  as 

(37) 

(38) 

(39) 

and  <r0  =  UM/w. 


^  eff  =  W  +  C±/D±> 


where 


r  ZUO 

c±  -  “T“ 


i  +  5 2  i  2r  ••  _  +  o-2 

1  0  2+p  71+4)  0  2+p(l+4) 


1+4  2  p  1  -  42  "l2  ,  ,2  f0  .  inT 

2  "o  2  2+p(l+4)  J  +6  [2  *  '^O  2 


D±=  1  -6  ±"-0 


In  most  cases  of  interest  the  loss  parameter  6  is  «  1  [compare 
Eq.  (15)].  For  o-Q  <  1  the  first  bracket  in  Eq.  (39)  does  not  vanish  anywhere 
in  the  interval  -1  <  4  <  I.  Under  these  conditions  C±  and  D±  can  be  replaced 
by  their  limit  for  6  -»  0  by  way  of  approximation. 


In  Fig.  5  the  imaginary  parv.s  of  the  two  circular  permeabilities  as 
calculated  on  this  basis  from  Eqs.  (37),  (38)  and  (39)  are  shown  as  functions 
of  the  relative  magnetization  for  =  0.  9  and  various  values  of  p.  The 

tensor  components  and  can  be  obtained  from  these  curves  in  accord¬ 
ance  with  Eq.  (12)  by  taking  half  the  sum  and  half  the  difference  of  p'|  and 
p"  e^.  Again  the  results  as  shown  are  based  upon  the  assumption  that  the 
core  domain  is  magnetized  down.  The  results  for  the  case  in  which  the  core 
domain  is  magnetized  up  are  obtainable  by  reversing  the  sign  of  4  and 
simultaneously  interchanging  p_,  and  p Comparison  of  Fig.  5  with  Figs.  3 
and  4  shows  that  the  imaginary  parts  of  the  permeability  show  more  hysteresi 
than  the  real  parts. 


It  is  interesting  to  note  that  the  absorption  coefficient  for  the 

negative  sense  of  circular  polarization  shows  a  peak  in  the  region  of 

positive  net  magnetization  (f,  >  0),  whereas  the  absorption  coefficient 

for  positive  polarization  shows  a  peak  in  the  region  of  negative  net 

magnetization  (£,  <  0).  Thus  the  absorption  in  general  peaks  for  the 

"anti-Larmor"  sense  of  circular  polarization,  i.  e. ,  that  sense  (defined 

with  respect  to  the  direction  of  the  net  magnetization)  which  does  not 

excite  the  ordinary  ferromagnetic,  resonance  of  a  completely  magnetized 

material.  This  result  is  in  agreement  with  the  experimental  observations 

2  5 

of  LeCraw  and  Spencer  and  Sandy.  Qualitatively  the  effect  can  be  attrib¬ 
uted  to  the  fact  that  a  small  down-domain  surrounded  by  a  large  up-domain 
has  a  high  resonant  frequency  (approaching  u^)  and  that  this  resonance 
is  excited  by  that  sense  of  circular  polarization  which  couples  best  t.o  the 
down  domain  (i.  e. ,  the  negative  sense  of  polarization  in  our  terminology). 

2 

This  is  substantially  the  qualitative  explanation  given  by  LeCraw  and  Spencer, 
which  is  confirmed  by  the  more  rigorous  theory  presented  in  the  present 
paper. 


In  Fig.  6  the  absorption  coefficients  for  the  two  senses  of  circular 

polarization  are  plotted  versus  t,  for  =  0.  5.  It  may  be  seen  that, 

compared  to  che  case  ioM/u)  =  0.  9  (see  Fig.  5),  the  absorption  is  now  much 

smaller.  1 1  particular  the  difference  between  and  p"_  (i.  e. ,  2k")  is- 

minute  compared  to  the  sum  (i.  e. ,  2p").  This  is  again  in  qualitative  agree - 

5 

ment  with  many  experimental  observations. 

For  a-g  >  1  it  is  no  longer  permissible  to  approximate  D±  by  its 
limit  for  <5  0,  even  when  6  «  1,  because  the  first  bracket  in  Eq.  (39) 

can  now  vanish  for  certain  values  of  t,  and  p.  This  is  the  Polder-Sinit 
resonance  mentioned  in  the  introduction  which  accounts  for  the  onset  of 
"low  field  loss". 

Some  examples  of  the  dependence  of  the  absorption  coefficients  upon 
the  reduced  magnetization  to  be  expected  under  these  conditions  are  shown 
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ion  coefficients  for  circular  polarization  ^  vs  reduced  mag- 
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positive,  broken  lines  to  negative  circular  polarization 


in  Figs.  7  and  8.  In  these  examples  p  has  been  chosen  as  4  and  6  as 
-2 

10  .  eff/6  versus  £,  is  shown  in  Fig.  7  for  o-q  =  0.  9,  1,  1.1  and  1.2. 

It  may  be  seen  that  near  o-q  =  1  the  absorption  increases  very  rapidly  with 
<Tq  for  almost  all  values  of  £,.  The  same  is  true,  though  to  a  lesser  extent, 
for  p'^  as  shown  in  Fig.  8. 

3c.  Application  to  Samples  Containing  an  Even  Number  of  Domains, 
Domains  of  the  Same  Type  Having  Equal  Volume 

Let  p  and  q  be  the  fractional  area  occupied  by  up  and  down  domains, 
respectively.  Thus 

P  +  q  =  1  (40) 

and  the  relative  net  magnetization  is 

i  =  p-q  •  (4i) 

It  is  shown  in  Appendix  2  that  for  a  cylindrically  symmetric  arrangement 
of  up  and  down  domains  the  components  of  the  effective  permeability  tensor 
can  be  expressed  as 


^eff  ,  2 

-  =1-0* 


2  y(n-) 


[ y(t r)  +  1]  2  -  >T2 


where 


k  ,  2  .  .2 

eff  1  -  it  -  yin-) 

-  =  cr - 

11  [yH  +  i] 


y(o')  =  q  fI/o‘)/pgN(o-) 


(42) 


(43) 


and  fj^(^)  and  are  defined  by  Eq.  (A  13)  in  terms  of  the  (1,1)  components 

of  matrices  S(tr)  and  T(<r)  given  by  Eq.  (A9).  For  N  =  2  and  4  these  matrices 


Fig.  7  Absorption  coefficient  for  negative  circular  polarization  p."  eff  vs  re¬ 
duced  magnetization  The  curves  apply  to  a  three-domain  configura 
tion  for  p  =  4,  6  -  10"  2,  -  0.  9  -  1.2,  and  Ha  =  0. 
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can  readily  be  evaluated  by  hand.  For  N  =  6,  10,  ?.0,  40,  80  they  were 
evaluated  numerically  on  a  computer.  For  simplicity  o r  was  taken 
as  real.  Thus  the  resultant  pg^  and  are  also  real.  The  results 
of  the  calculation  for  cr  =  0.  9  are  shown  in  Figs.  9  and  10.  The  case 
N  =  2  is,  of  course,  identical  with  the  three-domain  case  for  p  =  0. 
Figure  9  shows  the  dependence  of  p'e^  upon  t,.  With  increasing  N 
the  curve  becomes  more  nearly  symmetric  in  t,.  As  expected  from 
the  general  discussion  in  Section  3a  the  minimum  value  of  pe„  is 
independent  of  N. 

The  variation  of  ic'  ^  with  t,  is  shown  in  Fig.  10.  Also  shown 
in  this  graph  is  the  prediction  of  Rado's  theory.  It  appears  that  K'e^  as 
calculated  on  the  basis  of  the  present  theory  approaches  Rado's  value 
in  the  limit  of  large  N. 

The  remarks  concerning  hysteresis  of  the  variation  of  pe^  and 
Kejj.  with  f,  made  in  connection  with  the  three-domain  case  apply  equally 
well  in  the  present  case.  The  curves  shown  are  based  upon  the  assump¬ 
tion  that  the  core  domain  is  magnetized  down. 


4.  DISCUSSION 


The  theory  described  in  the  present  paper  differs  from  those  pre¬ 
viously  proposed  primarily  in  that  attention  is  first  focussed  on  the  deter¬ 
mination  of  the  local  rf  magnetic  field.  This  field  is  calculated  substantially 
rigorously,  subject  only  to  the  quasi-static  approximation.  In  order  to  carry 
out  the  calculation  it  was  necessary  to  assume  a  simple  domain  configuration. 
In  particular  the  explicit  calculations  of  the  effective  permeability  described 
in  Sec.  3b  and  c  are  based  upon  the  assumption  of  very  idealized  configura¬ 
tions.  Thus  the  results  of  Ihe  calculation  should  be  interpreted  with  appro¬ 
priate  caution. 

In  actuality  the  sample  will  contain  very  many  domains  of  different 
size,  and  they  may  be  oriented  in  many  different  ways.  For  such  a  domain 
configuration  the  resonant  frequencies  of  the  proper  modes  of  vibration  are 
very  likely  distributed  over  a  fairly  wide  frequency  band.  Thus  the  absorp¬ 
tion  coefficients,  considered  as  functions  of  magnetization  (£),  will  wery 
likely  not  show  sharp  peaks  as  deduced  on  the  basis  of  the  three-domain 
(see  Figs  7  and  8).  The  hysteresis  in  the  vs  t,  and  vs  £  curves 
may  be  expected  to  be  considerably  smaller  in  a  realistic  domain  config¬ 
uration  than  in  the  three-domain  configuration  assumed  in  Sec.  3b.  This 
is  confirmed  by  the  results  obtained  for  the  many  domain  model  discussed 
in  Sec.  3c  (compare  Figs.  9  and  10). 

For  a  comparison  with  experimental  data  some  allowance  must 
be  made  for  the  fact  that  the  domains  may  be  oriented  along  directions  other 
than  the  z-axis  as  assumed  throughout  the  present  paper.  The  general 
method  developed  in  the  present  paper  can  also  be  applied  to  such  domain 
configurations,  but  the  detailed  calculation  becomes  considerably  more 
involved. 


For  the  completely  demagnetized  state  the  random  orientation  of 
the  domains  can  easily  be  taken  into  account  in  an  approximate  fashion.  In 


this  case  the  off-diagonal  components  of  the  permeability  tensor  vanish. 

If  only  up-  and  down- domains  are  present  the  diagonal  components  cor¬ 
responding  to  the  x-  and  y-directions  are  given  by  Eq.  (30)  and  that  cor¬ 
responding  to  the  z-direction  equals  unity.  The  average  permeability  for 
random  orientation  may  then  be  expected  to  be  given  approximately  by  the 
average  of  the  permeabilities  applicable  in  the  three  principal  directions. 
This  yields,  according  to  Eq.  (30),  as  the  average  (isotropic)  permeability 
in  the  demagnetized  state 


^demag  3 


(  u/\) 


)2  -  (Ha  t^Mq)2-,!  j 

(Wy)2  -  H 2  -  2 


For  \(Ha  +  4ttM)/u  «  1  in  particular 


^demag 


1  UM  (tJM  +  2Ha) 


This  result  resembles  the  interpolation  formula  proposed  by  Bady,  which 
has  been  found  to  agree  quite  well  with  the  observations. 

The  permeability  defined  by  Eq.  (44)  is  real  for  u  >  y(4ttM  +  H  ) 

U  cl 

and  imaginary  for  yH  <  u  <  y  (4rMn  l-H. ).  It  thus  accurately  describes  the 
onset,  of  low-field  loss.  With  suitable  correction  to  allow  for  a  finite  line- 
width  and  possible  variations  of  the  anisotropy  field  from  one  domain  to  the 
other,  Eq.  (44)  has  been  found  to  agree  quite  well  with  the  observed  be¬ 
havior  of  (i1  and  1  in  the  demagnetized  state  for  frequency  between  approxi¬ 
mately  0.  3  GHz  and  10  GHz.  For  lower  frequencies  the  susceptibility  arising 
from  domain-wall  displacement,  which  is  not  taken  into  account  in  the  present 
theory,  apparently  becomes  important. 


The  theory  described  in  the  present  paper  accounts  qualitatively 
for  many  of  the  phenomena  observed  in  partially  magnetized  ferrites  at 
microwave  frequencies.  These  phenomena  include: 
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1.  The  onset  of  low  field  loss  near  u^/u  =  1, 

2.  The  variation  of  and  p^  with  net  magnetization,  in  par¬ 
ticular  the  hysteresis  observed  in  this  dependence, 

3.  The  variation  of  p^  and  p^  in  the  demagnetized  state  with 
uM/w  (see  Eq.  44), 

4.  The  variation  of  pT  and  p"  ^  with  net  magnetization,  in 
particular  the  fact  that  the  absorption  peak  occurs  on  the  anti- 
Larmor  side. 

5.  The  fact  that  for  £  0.  5  the  two  absorption  coefficients 

for  circular  polarization  become  nearly  equal,  i.  e.  , 

k  "  «  p"  ... 

eff  eff 

A  detailed  comparison  of  experimental  data  with  the  theory  will  be  presented 
as  a  separate  publication. 
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APPENDIX  II-A,  Derivation  of  Eq.  (2'9) 


Let  the  cross-sectional  area  of  the  n 


th 


domain  be 


a 


n 


(Al) 


and  assume  that  the  total  number  of  domains  (N)  is  even.  Then  according 
to  Eqs.  (24)  and  (28) 
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+  eff  _  1  al  Ai  +  a3  A3  +  -  •  •  +  aN_i  AN_i  "  [a2  A2  +  a4A4  +  "  '  aNAN^ 
^  a,  A,  +  a,  A.  +  -  •  •  +  a*T  n  A*T  .  +  [  a„  An  +  a,  A '  +•  •  •  a,T  A,  ' 
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N-l  N- 


4  4 


N  N 


Using  Eq.  (26)  this  result  may  alternatively  be  written  as 


(A  2) 


eff 

V- 


Sn(tr)  -  Tn(v) 
+  Tll(cr) 


where  S^tr)  and  T^cr)  are  the  (1,1)  components  (i.e. 
of  the  matrices 


(A3) 

first  row,  first  column) 


S(cr)  -  aj  l  +  a^  +•  •  •  +  a^_  ^  '  '  Si!?! 

T(tr)  =  [a2l  +  a4M3M2+’ •  •  +  aN  • -JVI2]  Mj  .  (A4) 


Using  the  definition  (27)  of  the  matrices  it  can  readily  be  shown  that 


S  n(<r)  =  F(cr) 

Tn(o ■)  =  (1  +  or)  G(cr)  (A5) 


where  F(cr)  and  G(cr)  are  both  even  functions  of  cr.  Thus  according  to  Eq.  (A3) 


APPENDIX  II-B,  Derivation  of  Eg.  (42) 


With 


aI  r  a3  =  ‘  ‘  '  "  aN-l  "  q,Tr2 


a2  a4  '  '  '  =  aN  ~  P  11  r2 


Equation  (A4)  can  be  written  as 


SM  =  qitr2  H  +  P2  +  P4P2  +•••  +  Pn.2'"^1 
T(r)  =  p»r^  [1  +P3  +  P5E3  +'"  +  Pn.1'--P31  M, 


where 


P  =  M  M  , 
~n  -  n  ~"n-  1 


Using  the  definition  (26)  of  one  obtains  for  integer  n 


-  n2  E 


1  -  cr 


(1  -  tr)o 


r2  n  (n  -  p) 


(1  +  or)  o'  r^  p  I  +  cr 


2  P 
n  -  p 


I  -  cr 


2  q 

n  +  q 


-  (1  +  a)  cr 


r2  n(n  +  q) 


^2n+l 


-  (1  -  <r)  o-  r2  q 


i  +  „2  a 
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(A  10) 


(All) 


If  we  define  f^(cr)  and  g^cr)  (symmetric  functions  of  <r)  by 
SnW  =  <5"r2  fNW 

Tn(<r)  =  prrr^  (1  +  tr)  gN((r)  ,  (A13) 

the  effective  circular  permeability  is  according  to  Eq.  (A3) 


where  y(o-)  is  defined  by  Eq.  (43).  Equation  (42)  is  now  obtained  by  using 
Eqs.  (11)  and  (12). 
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